Upper-critical dimension in a quantum impurity model: 
Critical theory of the asymmetric pseudogap Kondo problem 
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Impurity moments coupled to fermions with a pseudogap density of states display a quantum phase 
transition between a screened and a free moment phase upon variation of the Kondo coupling. We 
describe the universal theory of this transition for the experimentally relevant case of particle-hole 
asymmetry. The theory takes the form of a crossing between effective singlet and doublet levels, 
interacting with low-energy fermions. Depending on the pseudogap exponent, this interaction is 
either relevant or irrelevant under renormalization group transformations, establishing the existence 
of an upper-critical "dimension" in this impurity problem. Using perturbative renormalization group 
techniques we compute various critical properties and compare with numerical results. 



Non-trivial fixed points and associated phase tran- 
sitions in quantum impurity problems have attracted 
considerable interest in recent years, with applications 
for impurities in correlated bulk systems, in transport 
through nanostructures, and for strongly correlated lat- 
tice models in the framework of dynamical mean-field 
theory. Many of those impurity phase transitions occur 
in variations of the well-known Kondo model [1] which 
describes the screening of localized magnetic moments 
by metallic conduction electrons. A paradigmatic exam- 
ple of an intermediate-coupling impurity fixed point can 
be found in the two-channel Kondo effect. 

Non-metallic hosts, where the fermionic bath density 
of states (DOS) vanishes at the Fermi level, offer a differ- 
ent route to unconventional impurity physics. Of partic- 
ular interest is the Kondo effect in so-called pseudogap 
systems [2-8], described by the Hamiltonian 



H 



k 



Ja'S ■ s(0) 



(1) 



where s(0) is the conduction electron spin at the impu- 
rity site, S is the spin-i impurity, and the rest of the 
notation is standard. Most importantly, the fermionic 
density of states follows a power law at low energies, 
p{lu) = J2]i^i'^ ~ ^k) oc iVoli^r > 0)- Such a behavior 
arises in semimetals, in certain zero-gap semiconductors, 
and in systems with long-range order where the order 
parameter has nodes at the Fermi surface, e.g., p- and d- 
wave superconductors (r = 2 and 1). Indeed, in d-wave 
high-Tc superconductors non-trivial Kondo-like behavior 
has been observed associated with the magnetic moments 
induced by Zn impurities [9]. Note that the limit r — > oo 
corresponds to a system with a hard gap. 

The pseudogap Kondo problem has attracted substan- 
tial attention during the last decade. A number of studies 
[2-4] employed a slave-boson large- technique; signifi- 
cant progress and insight came from numerical renormal- 
ization group (NRG) calculations [5-7] and the local mo- 
ment approach [8] . It was found that a zero-temperature 
phase transition occurs at a critical Kondo coupling, Jc, 



below which the impurity spin is unscreened even at low- 
est temperatures. Also, the behavior depends sensitively 
on the presence or absence of particle-hole (p-h) asymme- 
try, which can arise, e.g., from a band asymmetry at high 
energies or a potential scattering term at the impurity 
site. A comprehensive discussion of possible fixed points 
and their thermodynamic properties has been given in 
Ref. 6 based on the NRG approach. 

So far, analytical knowledge about the critical proper- 
ties of the pseudogap Kondo transition is limited. Pre- 
vious works employed a weak-coupling renormalization 
group (RG) method, based on an expansion in the di- 
mensionless Kondo coupling j — NqJk- It was found 
that an unstable RG fixed point exists at j = r, corre- 
sponding to a continuous phase transition between the 
free and screened moment phases [2]. Thus, the per- 
turbative computation of critical properties within this 
approach is restricted to small r. Interestingly, the NRG 
studies [6] showed that the fixed-point structure changes 
at r = r* « 0.375 and also at r = i , rendering the relevant 
case of r = 1 inaccessible from weak coupling. Numerical 
calculations [6, 7] also indicated that the critical fluctua- 
tions in the p-h asymmetric case change their character 
at r = l: whereas for r<l the exponents take non-trivial 
r-dependent values and obey hyperscaling, exponents are 
trivial for r > 1 and hyperscaling is violated. These find- 
ings suggest to identify r = 1 as upper-critical "dimen- 
sion" of the problem (whereas r = plays the role of the 
lower-critical "dimension"). However, the relevant local 
degrees of freedom were not known so far. 

This letter will present the low-energy theory for the 
quantum transition in the p-h asymmetric pseudogap 
Kondo model, suitable for calculating critical properties 
close to and above r — \. The ingredients are a crossing of 
many-particle singlet and doublet levels, and their inter- 
action with the low-energy conduction electrons. Using 
RG methods, we show that this interaction is relevant for 
r < 1 and irrelevant for r > 1 . For r < 1 , we calculate the 
RG flow of the coupling, and compute universal critical 
properties by renormalized perturbation theory, in a sys- 
tematic expansion in (1— r). For r> 1, bare perturbation 
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theory is sufficient, and non-universal terms appear, as 
expected for a theory above its upper-critical dimension. 

Low-energy field theory. It is useful to consider the 
Kondo problem with a hard-gap DOS [10]: in the pres- 
ence of p-h asymmetry it shows a first-order quan- 
tum transition, i.e., a level crossing, between a Kondo- 
screcned singlet and a spin-i doublet state. As wc can 
understand the pseudogap DOS of (1) as consisting of 
an asymmetric high-energy part and a (asymptotically) 
symmetric low-energy part, wc can obtain an effective 
theory for (1) by coupling the above mentioned three 
(many-body) impurity states, obtained by integrating 
out high-energy degrees of freedom from (1), to the re- 
maining low-energy part of the conduction electron spec- 
trum. To represent the three impurity states we intro- 
duce auxiliary operators for pseudoparticles bg (for the 
singlet), fcr with a =t, i (for the doublet), with the con- 
straint bj&s + /^/(7 = Q = 1. The form of the coupling to 
the conduction electrons is dictated by SU(2) symmetry. 

The effective theory to be discussed in the remainder 
of the paper takes the form: 

H = {so + X)ftU + Xblbs (2) 

.A 

+ 5o(/>,c,(0)+h.c.) + / dklkl'-kcl^Ck^ 

J -A 

where we have represented the bath by linearly dispers- 
ing chiral fcrmions Cka-, and ^^(0) = J dk\k\^Cka is the 
conduction electron operator at the impurity site. The 
spectral density of the Ca{0) fermions follows the power 
law \uj\^ below the ultraviolet (UV) cutoff A. 

The parameter sq drives the system through the phase 
transition at sqc- For sq < soc the ground state is a dou- 
blet corresponding to a imscrccned local moment (LM), 
whereas so > soc leads to a screened singlet state (ASC, 
see below). For coupling go = the transition is simply 
a level crossing, and sqc = 0. As the original Kondo cou- 
pling Jk in (1) determines the relative positions of the 
singlet and doublet levels, we have (sq — sqc) oc {Jk — Jc) 
near the transition. 

The constraint for the Hilbert space of the pseudopar- 
ticles, (5 = 1, will be implemented using the chemical 
potential A ^ oo, such that obscrvablcs (O) have to be 
calculated according to (O) = \imx^oo{QO)x/{Q)x [11]. 

A crucial ingredient is the p-h asymmetry of the orig- 
inal model (1). It is clear that upon integrating out 
the high-energy part of the bath two many-body singlet 
states arise: in the strong-coupling limit they can be un- 
derstood as cither an electron or a hole bound to the 
impurity spin. Due to the p-h asymmetry of the under- 
lying model these two singlet states will have very dif- 
ferent energies, such that we can discard the high-energy 
state in the low-energy theory (2) [without loss of gener- 
ality we have kept the "electron" state] . The two singlet 
states are related by a global p-h transformation, and 
they constitute the ground states at the two asymmetric 
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FIG. 1: Schematic RG flow diagram for the effective theory 
(2). The horizontal axis denotes the energy difference be- 
tween doublet and singlet impurity levels, s, the vertical axis 
is the fermionic coupling g. (In the language of the infinite-?/ 
Anderson model the /-electron energy is e/ s s, and the hy- 
bridization V = g.) The thick lines correspond to continuous 
boundary phase transitions; the open circles are the critical 
fixed points, a) r > 1: g is irrelevant, and the transition is 
a level crossing with perturbative corrections, b) r < 1: g 
is relevant, and the transition is controlled by an interacting 
fixed point at g*^ = (1 - r)/3 + 0(1 - rf (4). 



strong-coupling fixed points (ASC) [6] of the pseudogap 
Kondo model. 

Interestingly, the theory (2) is formally identical to an 
infinite-[/ Anderson impurity model, and we will com- 
ment on this connection below. 

Scaling analysis and renormalization group. Tree level 
scaling analysis shows that dim[go] = = There- 
fore, r = l plays the role of an upper-critical "dimension" 
where go is marginal. The case r < 1 requires a RG 
treatment in analogy to the standard e expansion of tlic! 
(j)'^ theory. We shall employ here the familiar momen- 
tum shell method with cutoff A. For convenience, we 
introduce masses in each of the /(j and bg propagators, 
ruoa and TOos, keeping in mind that only their difference 
So = — TUos has physical significance. Dimension- 
less quantities arc defined as go = A*";^, moo- = Ania-, 
mos = ArUg. We set up RG equations for g, and rUg, 
where the fixed-point value of the non-linear coupling g"^ 
will be of order f . Subsequently, we can compute critical 
properties in a double expansion in g and f. 

To one-loop order we obtain the following RG beta 
functions: 

3 , 

P{9) = -rg+T;9\ 



/3(s) = -s + 3gh-g'' 



(3) 



where s = — nig. For r < 1 the trivial fixed point 

g* = s* = is unstable, and the critical properties are 
instead controlled by an interacting fixed point at 

2r ^ _ _2r 



9 = ^ : 



(4) 



At this fixed point, we find anomalous field dimensions 
rjg = 2g*^ , TTja = g*^ , and a correlation length exponent 



vz 



(l-3<?*')-i = l/r + 0(f2). 



(5) 



describing the vanishing of the characteristic crossover 
temperature in the vicinity of the critical point, T* oc 

|so - Sod"' [12]. 

The rcsuhing RG flow diagram is shown in Fig. 1. 
(No qualitative changes occur to two- loop order.) We 
have also indicated the stable fixed points reached for 
,s ±Qo: Large positive s implies a non-dcgcncratc sin- 
glet ground state of the impurity (ASC), where the inter- 
action go is clearly irrelevant. Large negative s describes 
a doublet impurity state (LM); for small go one can per- 
form the usual Schrieffer- Wolff transformation to project 
out the singlet state. This results in a pseudogap Kondo 
model of the form (1), with Jk oc 5o/I'^o| ~ here small 
Jk is irrelevant for r > [2, 13]. 

The critical fixed point (Fig. lb) shifts to larger values 
of g*^ , \s*\ with decreasing r. In the light of the numeri- 
cal results of Ref. 6 we expect that Eq. (2) together with 
an expansion in p, f yields the correct description of the 
critical properties for 0.375 « r* < r < 1. Using NRG 
we have numerically confirmed this expectation, i.e., the 
properties of the critical fixed point of the model (2) vary 
continuously as function of r for r* < r < 1. Considering 
that Eq. (2) represents an infinite-[/ Anderson model, 
we also conclude that in this r range the phase transition 
of the pseudogap Anderson model [6] is in the same uni- 
versality class as the one of the pseudogap Kondo model. 

We proceed by calculating a few observables in the 
vicinity of the critical point for r < 1 using stan- 
dard renormalized perturbation theory; diagrams are dis- 
played in Fig. 2. Power counting shows that the integrals 
appearing in universal quantities will display no UV sin- 
gularities, thus the cutoff A can be sent to infinity. The 
structure of the critical theory implies hyperscaling prop- 
erties, including lu/T scaling in dynamical quantities [12]. 
The local dynamic susceptibility follows a scaling form 



CO 



-Vx 
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T So - soc 



(6) 



where <&i is a universal crossover function (for fixed r), 
and Bi is a non- universal prefactor. A similar scaling 
forms holds for the local static susceptibility xioc (T) , and 
at the critical point Xioc{T) oc T^^+'^x. The anomalous 
exponent r]-^ can be calculated in an expansion in f, 



(7) 



which fits the numerical findings, e.g., 1 — r/^ — 0.928 ± 
0.002 for r = 0.9 [7]. Using hyperscaling, the results (5) 
and (7) are sufficient to determine all critical exponents 
associated with a local magnetic field [7] , in an expansion 
in f. In particular, the T ^ local susceptibility away 
from criticality obeys: 

Xloc(sO>SOc) oc {so - Socy , J = l^Z {1 - T]^) , 

Txioc{so<soc) oc (soc-so)'' , l' = vzr]^. (8) 
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FIG. 2: Feynman diagrams occurring in the perturbation the- 
ory for (2). Full/wiggly /dashed lines denote fa/hs/Ca prop- 
agators, a) Bare interaction vertex go- b) self-energy, c) 
ha self-energy, d) Diagrams entering the local susceptibility, 
needed to obtain the anomalous exponent rj-^ (7). 

As proposed earlier [6, 7] , limT->o Txioc can serve as an 
order parameter, as it vanishes continuously as Jk — *■ 
J~, and is zero for Jk > Jc- 

A scaling form similar to (6) holds for the conduction 
electron T matrix, describing the scattering of the c elec- 
trons off the impurity. At the critical point the T matrix 
will obey T{u:) oc 0;"^+''^. As in Ref. 13 we are able 
to determine an exact result for the anomalous exponent 
r]'rp [14], valid to all orders in perturbation theory. The 
argument is based on the diagrammatic structure of the 
T matrix: the relevant diagrams can be completely con- 
structed from full g vertices and full f /b propagators [13]. 
We obtain the exact result: 



Vt 



2f 



(9) 



Note that the similar result of Ref. 13 was derived for 
the p-h symmetric critical fixed point of the model (1), 
and is valid for < r < i. In contrast, Eq. (9) applies 
to the p-h asymmetric fixed point (4), and is valid for 
r* < r < 1. NRG calculations have found precisely this 
critical divergence, ImT(ti)) oc jt^'l"'", for both the sym- 
metric and asymmetric critical points [5, 15]. Notably, 
the T matrix can be directly observed in experiments, 
due to recent advances in low-temperature scanning tun- 
neling microscopy, as has been demonstrated, e.g., with 
high- temperature superconductors [15-17]. 

Lastly, we mention that it is also possible to compute 
the impurity residual entropy at the critical point, by 
expanding the free energy in g. The result is 



^imp(T = 0) =ln3 
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f+0(f2), (10) 



again in agreement with numerics. Fig. 14 of Ref. 6. 

Perturbation theory for r > 1. For bath exponents 
r > 1 the coupling go in the theory (2) is irrelevant in 
the RG sense, and therefore observables can be obtained 
by straightforward perturbation theory. Note that the 
cutoff A has to be kept explicitly, as integrals will be UV 
divergent. This also implies that hyperscaling is violated, 
and no u)/T scaling in dynamics occurs, as usual for a 
theory above the upper-critical dimension. 

We shall demonstrate the calculation of the local impu- 
rity susceptibility xioc- To lowest non-trivial order, xioc 
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is given by the convolution of two /o- propagators, cal- 
culated with the self-energy to second order in (/q. Note 
that no vertex corrections occur to this order due to the 
structure of the interaction. The self-energy, Fig. 2b, is 



l+r 



{2n) 
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-iu'-X-mos 



The real part of this self-energy is oc A and will renormal- 

ize the mass m,Qa-, this shift can be absorbed by defining 
ifia- = mo<T + R-eI]/^(A -|- fha-), similarly for rhg. Subse- 
quently, all quantities are expressed in terms of m^, fhs, 
and s = So — sqc = — fhs measures the distance to 
the transition point. 

At the transition, s = 0, the self-energy shows threshold 
behavior at T = 0, — Im5]/^(w -I- iri)/Tr cx gQCj''0(tj) with 
u! = u — X — fha- The low-energy behavior of the 
propagator follows as 

-ImG'/„(w + zr?)/7r = (1 - A)5{l;j) + B\CjY'-^<d{ui) 

with A,B oc 5q. The bg propagator has a similar form. 
When calculating xioc the T — > limit has to be taken 
with care, as the exponentially small tail in ImG/^ (w) at 
a; < contributes to xioc- One obtains for the imaginary 
part x(oc for T ^ and 1< r < 2: 



x!oc('^)A 



I - 2A S{aj)LO B 



+ - |a;r-^sgn(a;) . (11) 



The static local suceptibility at criticality is easily seen to 
follow Xioc{T) oc 1/T, implying rj^ = 0, in contrast to the 
result for r < 1 (7). The impurity residual entropy re- 
ceives no corrections from the coupling to the conduction 
electrons, thus 5i,np(T = 0) = ln3. 

Moving away from criticality, one of the fa-, bs propa- 
gators looses its 5{u)) contribution - for s > {Jk > Jc) 
this is Gf^, indicating that free- moment behavior is ab- 
sent in this regime. Consequently, the zero-temperature 
static local susceptibility is finite, but diverges upon ap- 
proaching the critical point according to xioc oc s^~'^. 
For s < 0, Xioc oc 1/T. Thus, the order parameter Txioc 
jumps at the transition point for r > 1. 

At r = 1, relevant for the d-wavc siipcrconductor, 
the quoted power laws remain valid, but are supple- 
mented by logarithmic corrections [3, 6, 15], as the cou- 
pling go is marginal in this case. Wc have found the 
leading corrections to be xioc oc l/(u; In^^"^ |w/A|) and 
ImT(a;)ocl/(|a;|ln^|a;/A|). 

Conclusions. We have developed an effective low- 
energy theory describing the quantum phase transition 
in the asymmetric pseudogap Kondo model - this was 
made possible by identifying the local degrees of free- 
dom relevant for the transition near its upper-critical 
"dimension". For bath exponents r < 1 we have ana- 
lyzed the interacting field theory using perturbative RG 



techniques, and established the existence of a transition 
obeying universal local non-Fermi liquid behavior and 
strong hyperscaling properties including w/T scaling of 
dynamical quantities. For r > 1 wc find a level crossing 
with perturbative corrections. Notably, the RG near the 
upper-critical dimension (r = l) needed to be formulated 
in variables completely different from the ones appropri- 
ate near the lower-critical dimension (r = 0). Our results 
are in excellent agreement with numerical data; they are 
applicable to impurity moments in unconventional su- 
perconductors [9, 15, 17] and other pseudogap systems. 
Further, the pseudogap Kondo problem permits a sharp 
transition upon application of a local magnetic field, this 
will be detailed elsewhere. 

We expect that field theories similar to ours can be con- 
structed for other impurity transitions, and will also be 
useful for the study of lattice models in dynamical mean- 
field theory and its extensions, with applications to lo- 
cal quantum criticality, discussed both in heavy-fermion 
metals [18] and high- Tc superconductors [19]. 

We thank R. Bulla, K. Ingcrscnt, M. Kircan, A. Rosch, 
S. Sachdev, Q. Si, and P. Wolfie for discussions. This 
research was supported by the DFG through the Center 
for Functional Nanostructures Karlsruhe. 



[1] 

[2] 

[s: 

[4] 
[5] 

[6; 

[7 

[s; 

[9 

[lo; 
[11 

[12 
[13 
[14 

[15 
[16 

[17 

[is: 



A. C. Hewson, The Kondo Problem to Heavy Fermions, 
Cambridge University Press, Cambridge (1997). 
D. Withoff and E. Pradkin, Phys. Rev. Lett. 64, 1835 

(1990). 

C. R. Cassanello and E. Fradkin, Phys. Rev. B 53, 15079 
(1996) and 56, 11246 (1997). 
A. Polkovnikov, Phys. Rev. B 65, 064503 (2002). 
R. Bulla et al., J. Phys.: Condens. Matter 9, 10463 (1997) 
and ibid 12, 4899 (2000). 

C. Gonzalcz-Buxton and K. Ingersent, Phys. Rev. B 57, 
14254 (1998). 

K. Ingersent and Q. Si, Phys. Rev. Lett. 89, 076403 

(2002) . 

D. E. Logan and M. T. Glossop, J. Phys.: Condens. Mat- 
ter 12, 985 (2000). 

J. Bobroff et al, Phys. Rev. Lett. 86, 4116 (2001). 
K. Chen and C. Jayaprakash, Phys. Rev. B 57, 5225 
(1998). 

P. Coleman, Phys. Rev. B 29, 3035 (1984). 
S. Sachdev, Quantum Phase Transitions, Cambridge Uni- 
versity Press, Cambridge (1999). 

M. Vojta and M. Kircan, Phys. Rev. Lett. 90, 157203 

(2003) . 

In the weak-coupling approach of Rcf. 13 the anomalous 
exponent of the T matrix is defined differently. 
M. Vojta and R. Bulla, Phys. Rev. B 65, 014511 (2002). 

E. W. Hudson et al, Science 285, 88 (1999); S. H. Pan 
et al. Nature 403, 746 (2000). 

A. Polkovnikov et al, Phys. Rev. Lett. 86, 296 (2001). 
Q. Si, S. Rabello, K. Ingcrscnt, and J. L. Smith, Nature 
413, 804 (2001) and Phys. Rev. B 68, 115103 (2003). 
[19] C. Panagopoulos et al, Phys. Rev. B 66, 064501 (2002). 



